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ABSTRACT 

Two pressure chambers, for compression experiments with liquids from zero 
to 2.2 GPa pressure, are described. The experimentally measured compressions 
are then compared to theoretical values given by an Isothermal model of 
equation of state recently Introduced for solids. The model describes the 
pressure and bulk modulus as a function of compression for different types of 
lubricants with a very high accuracy up to the pressure limit of the high 
Pressure chamber used (2.2 GPa). 

In addition the Influence of temperature on static solidification pressure 
was found to be a simple function of the thermal expansion of the fluid. 
NOMENCLATURE 

2 

B bulk modulus, N/m 

2 

B q bulk modulus at zero pressure, N/m 

m mass, .kg 

2 

P pressure, N/m 

t temperature, °C 

. 3 

v volume, m 

3 

v .j volume at the solidification pressure, m 
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v sol 0 ''°^ ume *be lubricant should become If It were decompressed from 

the solid state down to zero pressure keeping the solid state 
properties, m^ 

x parameter In bulk modulus and pressure equation 

x sol (v/v sol,0 )1/3 

a coefficient of thermal expansion, °C _1 

y dT^/dp, shear strength Increase with pressure 

6 coefficient of thermal volume expansion, °C ' 

3 

p density of lubricant, kg/m 

n parameter In bulk modulus and pressure equation 

2 

t l shear strength of lubricant, N/m 

INTRODUCTION 

In the calculation of the oil film thickness In elastohydrodynamlc 
lubrication, the compressibility and the change of lubricant behavior from 
liquid to solid type, are very Important. The compressibility of the 
lubricant determines the shape and the thickness of the oil film In the 
central almost parallel part of It. Lubricants with high compressibility give 
a thinner central oil film than stiff lubricants. If the compressibility Is 
high enough, the minimum film thickness will not appear at the outlet of the 
elastohydrodynamlc contact, but In the center of the contact. There the area 
of low film thickness will be much larger than the outlet area with low film 
thickness. Due to the large close contact area and the high pressure In this 
area, the heating and energy dissipation, If the oil film falls, Is much more 
severe than It should be If the minimum film thickness were the same, but only 
occupying the side lobes and outlet part of the contact. 

i 

The change of solidification pressure with temperature Is also Important 
for the lubricant film thickness but Is even more Important for the power 
dissipation In the elastohydrodynamlc contact. If the temperature rise Is 
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high enough that the oil behaves like a liquid, even In the central part of 
the contact, the friction forces will decrease and become lower than the 
limiting shear stress of the lubricant. In this paper two test apparatuses 
are described, one for measurement of the bulk modulus at low pressures and 
the other for measurement of the compression curves for high pressures up to 
2.2 GPa for liquid lubricants. A physical model, built upon consideration of 
the atomic Interaction forces In the lubricants, describes the 
pressure-compression relationship, for all lubricants tested, using only two 
parameters. These parameters are the bulk modulus at zero pressure B q and 
the pressure derivative of the bulk modulus at zero pressure (3B o /dp)p_g. 

The high pressure chamber was also used to measure the change In 
solidification pressure with a change In temperature. 

The High Pressure Chamber 

A detailed description of the high pressure chamber can be found In 
Ref. 1. Equipped as In Ref. 1, the high pressure chamber can be used to 
determine the solidification pressure for lubricants by measuring the shear 
strength of the fluid as a function of pressure. A description of the 
specially equipped high pressure chamber for high accuracy measurements of the 
compression of lubricants can be found In Ref. 2. 

A drawing of the apparatus Is shown In Fig. 1. As Is described In Ref. 2 
the pressurized volume between the plungers can be measured very accurately. 
The maximum error being about 2x10 of the total volume. It Is not only 
possible to measure the compression of the fluid In the high pressure chamber 
but also to measure the shear strength of the fluid as a function of the 
pressure. In these experiments the solidification pressure Is said to be the 
pressure at which the shear strength of the fluid starts to Increase linearly 
with pressure as a^/dp = y. 
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The high pressure chamber was used for measurement of the compression of 
the fluids from 0.4225 to 2.20 GPa. At pressures below 0.4225 GPa some of the 
lubricants with low viscosity at these pressures did leak past the steel seals 
as the steel seals were not plastically deformed to seal perfectly between the 
plungers and the Inner cylindrical surface of the high pressure chamber. 
Apparatus and Method to Determine 

Since leakage was observed at the lower pressures. It was not possible to 
obtain a reliable value of the bulk modulus at zero pressure (B q ). Therefore 
a new apparatus was built to measure B q . It consisted of a closed volume 
with very thick steel walls and a pressure transducer connected to It. 

The volume was then completely filled with the test lubricant and sealed. 

A temperature regulation system and a heating coll around the steel 
cylinder made It possible to keep the temperature and the pressure due to the 
thermal expansion of the oil at predetermined values. The density of the oil 
In the steel cylinder was 


m 


where m was the mass and v was the volume of the lubricant. 
The density change due to a temperature change was 


3p m 3v _ p 3v 
3t = ” v 2 3t * “ m 3t 

3v m 3p 

3t = ' 2 3t 

P 


( 1 ) 


( 2 ) 

(3) 


1 3V _ 1 d£ 

v at = " p at 

The relative Increase of the steel cavity volume caused by the temperature 
Increase was 3aAt where a = 11x10^ °C - ^ was the coefficient of thermal 
expansion for the steel. 


(4) 
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A temperature Increase of At, created a compression of the lubricant 
determined by the thermal expansion of the oil minus the thermal expansion of 


the container minus the elastic Increase of the container volume. 

The thermal expansion of the lubricant due to the temperature Increase 
was 8v/at»At and the thermal expansion of a unit volume of the lubricant 
was (3v/3t)/v»At. The thermal expansion of the unit volume of the steel 
container was 3av At/v = 3a At. The elastic volume Increase of the steel 
container due to a pressure Increase Ap was Av/Ap'Ap and the relative 
elastic volume Increase of the container was (l/v)»(Av/Ap)Ap. The 
compression of the lubricant was (-3v/3p)»Ap and the relative compression 
of the lubricant caused by the pressure Increase AP was (l/v)(-3v/3p) »Ap. 


(I % ' 3 “) 


1 AV 
V Ap 


or If Eq. (4) Is used 


1 _ 

B_ 


1 3V_ 1 AV 
V dp = V Ap 


1 3v 

7^ * Ap 

(5) 


(6) 

:/ Ap 

3p\ At 
3t/ Ap 

(7) 


For the apparatus (l/v)(Av/At) = 0.015 GPa - ^ and 3a = 33x10 °C _1 . 

For the different fluids the values measured at 20 °C are given In 
Table I. The values of p and 3p/3t were measured using a precision scale 
and calibrated volumes. 

Isothermal Volume Change of Lubricants With Pressure 

Vlnet, Ferrante, Smith, and Rose (3) calculated the pressure as function 
of compression for Ionic, metallic, covalent, and rare gas solids. Their 
equation of state gave the pressure 

3B 

p = -f- (1 - x) exp(n(l - x)) (8) 
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and the bulk modulus 


B * [2 + (n - l)x - nx 2 J exp(n(l - x)) 

x 

as a function of the relative compression x = (v/Vq) 1 / 3 
Here B 0 was the bulk modulus at zero pressure and 


(9) 



was a function of the pressure derivative of the bulk modulus at zero pressure. 
These equations are valid for solids at Isothermal conditions without phase 
change. 

As the Interatomic forces are similar In solids and liquids In 
compression, the equation of state for solids was applied as an approximation 
for the equation of state for liquid lubricants. If the lubricants behave In 
the same way as solids, the experimental points will fit the pressure 
deformation curve given by Eq. (8). 

Lubricant molecules are normally large and complex and have different 
types of bondings. This Is why eventually B q and n could vary during 
the compression. The lubricants tested were known to convert from liquid to 
solid behavior within the pressure range tested. As the pressure range for 
the test apparatuses described earlier did not overlap, the connection between 
low pressure and high pressure behavior was done using curve fit techniques. 
The pressure and bulk modulus equations were assumed to be applicable to 
lubricants and the constants In the pressure equation was determined by curve 
fitting of the measured data. 

Compression Through the Solidification Boundary 

For the liquid state the Eqs. (8) and (9) were used up to the 
solidification pressure. In the solid state the same types of equations were 
used, but the values of n, B q , and x were changed. 
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( 11 ) 


3B 

P = -p" 1 (1 - x $$ ) exp(n (1 - x ss )) 
x ss 

Here, B qs was the bulk modulus for the solid state at zero pressure. It 
could not be directly measured, but It was calculated using the equation 


B = -f-[ 2 + (n s - 1) x $s - n s xOexp (n s (l - x $s )) (12) 

x ss ^ J 

and the assumption that the bulk modulus was a continuous function of 
pressure, which means that B for the solid state and B for the liquid 
state at the solidification pressure was the same. At the solidification 
pressure B = B s » p = p s , and x $$ = x sq1 

36 

P S - -f (1 - x sol> exp <"s (1 - x sol» < 13 > 

x sol 

The pressure In the solidified state can be written as 


P = 


3(1 - x ss ) 


p s x sol 


x 2 exp(n s (l - x $s )) 3(1 - x sol ) exp(n s (l 


Sol” 


(14) 


g- • n‘ ‘ e * p( V»i (1 ' x «» )) 

s "os (1 - Sol* 


(15) 


1 /3 

where x^^ = (v/v ,) was the relative volume calculated from the sol 1 d 1 - 
os sol 

flcatlon volume raised to the power 1/3. 

x $o ^ was the relative volume at the solidification pressure raised to 
the power 1/3 assuming solid type behavior all the way down to zero pressure. 


x = x • x . 
ss os sol 


( 16 ) 


Curvefitting the experimental points gave x .j and n s 


Analytical Solution for x 


sol 


It Is possible to solve for x^ If B q , n^, n s , B $ , and p $ 
are known. The equations for the pressure as a function of relative volume 
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change are Eq. (8) for pressures up to the solidification pressure p $ and 
Eq. (15) for pressures above p $ . 

x SQ j can be solved for by using the fact that the bulk modulus Is the 
same at the solidification pressure for both the liquid and solid state. 

x sQ i Is then a function of the solidification pressure p s , the bulk 
modulus B $ , and the parameter n $ for the solid state. 

£2 + (n L - l)x - n L x 2 J exp(n L (l - x)) 

- * <*s - ’) ’■sol - "s x soll exp( "s (1 - x sol> (17) 

x sol - J 


But from Eq. (13) 



If ns = 0» as the experimental measurements suggest. 





( 21 ) 


( 22 ) 
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and 


36 os “sol 

"s = 1 - -SO! ‘ 

x $o ^ can be used to determine the temperature Influence on the 

solidification pressure, as will be shown at the end of this paper. 

Numerical Calculations of B and n. 

o 1 

In the equations for p and B the values of x and n are not known. 
The compression of the fluid Is only known relative to the volume at the 
Initial pressure 0.4225 GPa. The equation for p will contain this unknown 
compression. 

If the volume at zero pressure Is v q and the volume at the starting 
pressure Is v $ the equation for the liquid state pressure Is the same as 
Eq. (8) where x Is now given as 




>1/3 


yl/3 

= 'd x s 

5 > 


0/ vs o/ Vs/ VO. 
v/v s Is the relative volume calculated from the starting pressure 

\l/3 


(24) 


•fe 


) 


X - X X 

s OS 


(25) 


Using this expression for the relative volume, the pressure, and bulk modulus 
equations will be 


36 


P = ~ 2~2 (1 " x o x s ) ex P( n L (1 " x o x s» 


B = B. 


xx 

os 


t "l - 1 


2 2' x x ’ ” n L 
X o x s 0 S 


exp(n L (l - x Q x s )) 


(26) 


(27) 
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Here only the two values of x $ and n L are unknown. They can be solved 
for by least squares techniques using the measured values of p and x q . 
Results From Curveflttlnq of the Experimental Heaurements 

As the values of B q and n are different In the liquid and solid 
states for the lubricants. It Is possible to find the solidification pressure 
for a lubricant using curvefitting techniques. The solidification pressures 
found In this way, just by observing the change In curvature for the 
pressurevolume curve, are compared to the solidification pressures found by 
measuring the shear strength of the lubricants In Table II. The two 
completely different methods to determine the solidification pressures give a 
good agreement. The mean absolute difference being about 5 percent and the 
maximum difference being 8.9 percent. The solidification pressure In the 
curve fitting experiments has been assumed to coincide with one of the 
measured pressures. The pressure step from one experimental point to the next 
at the solidification pressure could be as much as 15 percent though, 
explaining most of the scatter. In Table II are also given the measured or 

3 

calculated parameter values for the different lubricants. v $ tart^ v o = x s 

Is the relative volume of the lubricant at the start pressure, 0.4225 GPa, of 

the experiment. t> l Is the slope parameter describing the slope of the 

curves In Fig. 3. The parameter n $ for the solid state of all tested 

lubricants was found to be zero, v ,/v . . Is the volume at the 

sol start 

solidification pressure divided by the volume at the pressure 0.4225 GPa. If 
the stiffness of the oil was the same as at the solidification pressure and 
the parameter n $ = 0 from zero pressure to the solidification pressure, 
the fictitious relative volume at the solidification pressure should be 
v so -|/v so i q. This relative volume, also given In Table II, will be used 
In the calculations of the change of solidification pressure with temperature. 
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In Fig. 2 the compression as a function of the pressure Is shown for six 

different fluids. The solid lines are the theoretical curves from Ref. 3 

fitted to the experimental points shown In the figures. The fits between the 

theoretical curves and the experimental points are remarkably good. In Fig. 3 

the natural logarithm of the nondlmenslonal bulk modulus Is shown as a function 

1 /3 

of the compression parameter x = (v/v Q ) . The theoretical solid lines 
and the experimental points are shown for the six fluids. The parameter n 
In Eq. (8) Is shown In this type of plot as the slope of the curves. It Is 
clearly visible from all six diagrams In Fig. 3 that the parameter n $ for 
the solid state Is zero. 

The Change of Solidification Pressure Caused by a Temperature Change 

As earlier mentioned Vlnet, Ferrante, Smith, and Rose (3) published 
equations for the pressure as a function of compression and the bulk modulus 
as a function of pressure. These equations were valid for Isothermal 
conditions without phase changes. Experiments with lubricating oils showed 
that the bulk modulus at low pressures Is almost Independent of the 
temperature. The total pressure and compression of an oil needed to convert 
It from a liquid to a solid state Increases with Increasing temperature. The 
Idea behind this part of the paper Is that the absolute density (or volume) 
determines the switch from liquid to solid state for lubricants. This means 
that the thermal expansion due to a temperature Increase has to be compensated 
by a pressure Increase which Increases the density to the value needed for 
solidification. 

If the volume coefficient of thermal expansion for the fluid Is 4 °C~\ 
the change In solidification pressure with temperature Is given by the thermal 
expansion 

vA • At = 3 Axv (28) 

where 
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V + AV 
V 


3 /, 

x = (1 


+ Ax) = 1 + 3 AX 


2 3 

+ 3Ax + Ax 


(29) 


~~ « 3 Ax since Ax « 1 (30) 

where x Is the parameter In the pressure equation. The pressure as a 
function of B q , x, and n Is given by Eq. (8). 

The experimental measurements In Fig. 3 show that the parameter value n 
In the solid state Is n $ = 0. This Indicates a simple method to calculate 
the change In solidification pressure caused by a change In temperature. 

At the solidification pressure for a lubricant with the temperature t^ 

Eq. (8) can be rewritten as 


In 


p sl • x sol 
3 < , -*sol> B os 


’ (1 - SolN ■ 0 


(31) 


3 

where x SQ ^ = v so i/ v so ] q given In Table II. 

If the lubricant then Is further compressed to the pressure p $2 and 
heated to the temperature t^, giving the thermal expansion 3 Ax = 6(t 2 - ) 

from Eq. (28), then Eq. (8) will give 


In 


s2 


< x sol - 4X) 


3(1 - x $ol * 4x)B 


OS 


= (1 - x so1 ♦ Ax)n s = 0 


(32) 


The solidification pressure at temperature t 2 divided by the solidification 
pressure at temperature t 1 can then be solved for from Eqs. (31) and (32). 

^2 , 


_s2 

J sl 


sol 


x sol + Ax 


x sol - AX - 


1 - X 


sol 


(33) 


Example 

For a typical lubricating oil n L Is of the order 10 to 30 and n $ « 0. 
The compression needed to solidify It Is Av/v q = 0.2 + 0.05 and If the oil 
had had the same stiffness as In the solidified state the compression needed 
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to reach the solidification pressure would have been (av/v 0 ) S0 ^ = 0.040 - 
0.084. The thermal expansion caused by a temperature Increase of 20 °C (from 
20 to 40 °C) Is from Table I. 

ax = 7.25X10" 4 x 20/3 = 0.00483 

3 

and from Table II x^ = 0.9359, x $Q ^ = 0.978 for the lubricant 
dltrldecyladlpate 26 cSt. 


This gives 






- tx > 2 

1 - "sol * 6X 



x so1 4 *» 
" x sol 


0.978 
10.978 - 0.00483^ 


0.02683 

0.022 


(34) 


(35) 


*51 


= 1.232 


(36) 


The compression experiments gave for the two temperatures 40 and 20 °C, 
P s2 /P sl = 1-95/1.58 « 1.234. 

For a 450 cst polyalphaolef Ine lubricant x $1 = 0.971, ax = 0.003937 for 
at = 20 °C gave p s2 /p sl = 1-14? and the compression experiments gave 
P s2 /p sl = 2 * 06/1 * 80 = 1-144. 

For a naphthenic base distillate with a very low solidification pressure 
(p $ = 0.67 GPa at 20 °C) the constants from the numerical curve fitting was 
x sq i = 0.9865 ax = 0.00896 for a temperature change of 40 °C. 



( 0.9865 V 
y0. 97754 J 


0.02246 

0.01350 


1 .694 


(37) 


and the experiments gave P^/Pji = 1-13/0.67 = 1.687. 

These values show a remarkably good agreement between theory and 
experiments. 
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CONCLUSIONS 


Two pressure chambers, for measurement of compressibility and bulk modulus 
for liquids from zero to 2.2 GPa pressure, are described. The experimentally 
measured compressions are then compared to values given by an Isothermal 
equation of state for solids In compression. The equation describes the 
pressure and bulk modulus for different types of lubricants with a very high 
accuracy up to the pressure limit of the high pressure chamber used (2.2 GPa). 

The Influence of temperature on static solidification pressure was found 
to be a simple function of the thermal expansion of the fluid. 
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TABLE I. - EXPERIMENTAL RESULTS 


Lubricant, 

nr 



p 

3p 

At 

1 3 V 


at 

A P 

v ap 

kg/m 3 

kg/m 3 °C 

°C/6Pa 

GPa -1 

846 

-0.591 

1043 

-0.679 

910 

-.660 

910 

-.615 

1005 

-.726 

995 

-.672 

1004 

-.716 

910 

-.604 

931 

-.626 

940 

-.586 

892 

-.629 

910 

-.597 



TABLE II. - MEASURED AND CALCULATED VALUES 



Solidification Solidification Bulk v s tart 

pressure from pressure from modulus — 

curve fit, shear strength at zero o 
GPa measurement, pres- o 

GPa sure B , x' 3 

o s 

GPa 


1.650 1.682 1.473 0.8654 

1.449 1.561 1.626 .8697 

0.995 1.092 1.488 .8736 

1.124 1.213 1.656 .8734 

.676 0.706 1.706 .9014 

.834 .839 1.675 .8945 



13.65 


13.47 


14.84 


13.71 


25.27 


20.93 


v sol 

< 

in 

O 

v start 

v sol ,0 


x 3 
sol 

0.8718 

0.9161 

.8899 

.9359 

.9252 

.9581 

.9247 

.9525 

.9759 

.9600 

.9564 

.9501 


Lubricant 

type 


fine 450 est 

Ditndecylade- 
pate 26 cst 

Polypropylene- 
glycole 175 cst 

3 olypropylene- 
glycole 80 cst 

Naphthenic base 
distillate 26 cst 

Naphthenic base 
raffinate 23 cst 
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